Superconducting gap structure of Celrlns from field-angle-resolved measurements 

of its specific heat 
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In order to identify tlie gap structure of Celrlnj, we measured field-angle-resolved specific heat C(0) by 
conically rotating the magnetic field H around the c axis at low temperatures down to 80 mK. We revealed that 
C(<f)) exhibits a fourfold angular oscillation, whose amplitude decreases monotonically by tilting H out of the 
at plane. Detailed microscopic calculations based on the quasiclassical Eilenberger equation confirm that the 
observed features are uniquely explained by assuming the (f,.2_,,2-wave gap. These results strongly indicate that 
Celrlnj is a d^.2_,,2-wave superconductor and suggest the universal pairing mechanism in CeMInj (M = Co, Rh, 
and Ir). 

PACS numbers: 74.70.Tx, 74.25.Bt, 74.25.0p, 74.20.Rp 



The heavy-fermion systems CeMIns (M - Co, Rh, and 
Ir) have been extensively studied because they exhibit uncon- 
ventional superconductivity near the antiferromagnetic (AF) 
quantum critical point (QCP). Especially, much effort has 
been expended to identify the superconducting (SC) gap struc- 
ture, a challenging issue that is closely related to the iden- 
tification of the pairing mechanism. Recently, field-angle- 
resolved experiments along with theoretical works have con- 
firmed that CeCoIns is a c/f2_,,2-wave superconductorji^ Be- 
cause of the t/j-2_,,2-wave gap and the proximity of the sys- 
tem to the AF state, it is now widely accepted that the pairing 
interaction in CeCoIns is AF spin fluctuations. Pressure in- 
duced superconductivity in CeRhlns has also been examined 
by field-angle-resolved specific heat measurements and it is 
argued to have the same pairing symmetry ii 

By contrast, the possibility of a different pairing mechanism 
has been suggested for Celrlnj due to its unusual behavior of 
the SC phasc/i. The transition temperature T^, which is 0.4 K at 
P=Q, increases to 0.8 K under a high pressure of P-2. 1 GPa^ 
although strong AF fluctuations existing at ambient pressure 
are rapidly suppressed by increasing P as the system is further 
pushed away from a hypothetical AF QCP. Moreover, when 
Rh is doped to Ir sites, shows a cusp-like minimum before 
reaching the maximum Tc of ~1 K around the onset of an 
AF statei^i^ These observations suggest the presence of two 
distinct SC phases in CeRhi- JrJns. 

While several attempts have been made to uncover the 
gap structure of Celrlns, two conflicting possibilities have re- 
mained. Kasahara et alX reported a fourfold angular oscil- 
lation in the thermal conductivity /c(0) when H is rotated in 
the basal plane, and attributed its origin to the vertical line 
nodes of the ii(.2_,,2-wave gap. On the other hand, Shakeripour 
et ali^^^ examined the effect of impurity scattering of k(T) 
for a current parallel and perpendicular to the c axis, and pro- 
posed the gap function of Celrlns to be either k, or k.{k^ -t- ik,,), 
both of which have a horizontal line node only on the equator, 
in sharp contrast to the c/j.2_,,2-wave gap. At present, neither 
of these two possibilities can be ruled out. The temperature 
variation of the anisotropy KdKa in Ref. [l^ can be explained 



within the dyp.^^p. symmetry if the Fermi surface has small devi- 
ations from the cylindrical symmetry^* whereas the fourfold 
oscillation in a-(0) can be explained by the horizontal line node 
gap if an in-plane anisotropy of the effective mass and/or min- 
ima exist in the azimuthal variation of the gap amplitudeJ^ 
Thus, the gap structure of Celrlns has remained controversial. 

In order to settle the controversy over the gap structure 
of Celrlns, we have performed field-angle-resolved specific 
heat C((^, 6) measurements down to 80 mK. The C(0) mea- 
surements, where (p denotes the in-plane azimuthal angle of 
H, have proven to be quite useful to determine the direction 
of nodes in the momentum space of bulk superconductors 
Here we extend the method to measure the polar angle 6 de- 
pendence of C(0) by which a detection of the horizontal line 
node can be made. We revealed that C(0, 6) exhibits a clear 
fourfold oscillation as a function of with H rotated around 
the c axis, and its amplitude is monotonically suppressed by 
tilting H out of the ab plane (decreasing from 90°). The 
results are compared with the microscopic theory by solving 
the quasiclassical Eilenberger equation self-consistently, and 
are found to be in good agreement with the iiv2_y2-wave gap, 
but are in a sharp contrast with the behavior predicted for the 
horizontal line node gap. 

The single crystal of Celrlns used in the present study 
(Tc = 0.4 K, 40.4 mg) was grown by the self-flux method. The 
specific heat was measured by the relaxation and the standard 
adiabatic heat-pulse methods in a dilution refrigerator (Ox- 
ford Kelvinox AST Minisorb). Magnetic fields were applied 
in the xz plane by using a vector magnet consisting of hori- 
zontal split-pair (5 T) and vertical solenoid (3 T) coils. By 
rotating the refrigerator around the z axis using a stepper mo- 
tor mounted at the top of a magnet Dewar, three dimensional 
control of the magnetic field direction is achieved. We con- 
firmed that the addenda contribution was always less than 5% 
of the sample specific heat and had no field-angle dependence. 
An accurate (< 0.1°) and precise (< 0.01°) field alignment 
with respect to the crystalline ab plane was accomplished 
by making use of the C{0) data which reflects the tetragonal 
anisotropy of Hc2- 
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FIG. 1: (Color online) Magnetic field dependence of the nuclear 
subtracted specific heat divided by temperature C^IT for (&) H \\ 
[100] and Qo) H \\ [001]. The solid lines represent fits to the data by 
a VH + b. The inset in (b) shows the temperature dependence of the 
upper critical field determined from the present study. 



The specific heat of Celrlns at zero field is known to exhibit 
an upturn on cooling below about O.IK due to a quadrupole 
splitting of the "^in (/ = 9/2) and 'yi i'^^jj. ^ ^ 3/2) nuclear 
spins.ii^^ In this paper, the nuclear Schottky contribution was 
subtracted from the data assuming Cn = (ao + a\H^)IT^, 
where ciq was adjusted so that the resulting electronic con- 
tribution Ce - (C — Cn) at // = becomes proportional to 
at low T . The value of oq (= 0.38 mJ K/mol) thus ob- 
tained was in good agreement with the calculated one (~ 0.34 
mJ K/mol) using the parameters determined from the "^In 
nuclear quadrupole resonance experiment, ' - The coefficient 
a\ was calculated from the nuclear Zeeman splitting. It can be 
shown that Cn is independent of the in-plane field orientationji 
Because of the nuclear Schottky contribution, the C{(f>, 6) mea- 
surements were limited to the temperature range T > 80 mK. 

Figures 1(a) and 1(b) show C^{H) for || a and H \\ c, 
respectively. As represented by the solid lines, Ce{H) is 
proportional to the square root of H at low H and low T. 
This behavior supports the presence of line nodes in the SC 
gap. The C^iH) behavior near H^i is in good agreement with 
the calculated result for a two-dimensional d-wave supercon- 
ductor with a relatively small Pauli paramagnetic parameter 
fi - 0.86,- At 80 mK for H \\ a, a cusp-like structure is 
observed at 0.06 T, which might originate from the multi- 
gap superconductivity, as reported in Sr2Ru04 (Ref. [l9h and 
MgB2^ 

Based on the Ce{H) data, we determined the upper criti- 
cal field Hc2{T), as plotted in the inset of Fig. 1(b). At 
the lowest temperature 80 mK, //c2 of the present sample is 
1.0 T for H II a and 0.55 T for H \\ c. The reduced ra- 
tio a - -Hc2/{TcdHi^2/<iT\T=Tj is estimated to be 0.47 for 
H \\ a and 0.54 for H \\ c, which agree well with the results 
of Ref. |2l[ These results also suggest that the Pauli paramag- 
netic effect in Celrlns is relatively weak compared with that 
in CeCoIns [a = 0.26 (Ref.IH]. 




FIG. 2: (Color online) Variations of C^/T as a function of the az- 
imuthal angle between the [100] axis and the magnetic field applied 
in the ab plane (6 = 90°). The dashed lines represent fits to the data 
by Co(T) + Ch{T, H)(l - A4 cos 40). 



Figure 2 shows the (p dependence of Cg/T obtained by ro- 
tating H in the ab plane, where (p is measured from the [100] 
direction. For each data point, the Ce(0) value was determined 
by an average of ten successive measurements. The error bar 
is estimated to be several mJ/mol K- for T > 1 10 mK, while 
it becomes about 10 mJ/mol-K^ at 80 mK due to the large nu- 
clear contribution. A clear fourfold oscillation was observed 
in Ce(0) in the wide T and H region. We carefully confirmed 
the absence of a twofold oscillation in Ce((p) which guaran- 
tees the accurate and precise alignment of the magnetic field 
with respect to the ab plane. In this intermediate T regime 
(0.2 < r/Tc < 0.5), Cg{(p) becomes minimum in fields along 
the {100) directions. 

To characterize the H and T variations of the fourfold os- 
cillation, we fit the data to the expression 

Ce(r, H, 0) = Cq(T) + Ch(T, H)(l - A4 COS40) (1) 

as represented in Fig. 2 by the dashed Unes. Here, Co and Ch 
are the zero-field and field-dependent components of Cg, re- 
spectively, and A4 is the amplitude of the fourfold oscillation 
normalized by Ch- Figure 3(a) shows the T dependence of A4 
at fioH^O.09, 0.2, and 0.5 T. We found that A4{T, H) exhibits a 
peak at around 0.3rc and 0. 15//c2, and rapidly decreases down 
to zero around 0.2rc. This is a feature that, to the best of our 
knowledge, has not been observed in previous /c(0) measure- 
ments done mainly above 0.4rc,— and strongly suggests the 
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FIG. 3: (Color online) Temperature dependence of the normalized 
fourfold amplitude A4 at several fields obtained by (a) the present ex- 
periment and (c) the microscopical calculations assuming the d^i^yi- 
wave gapi^ (b), (d) Contour plots of A4(T, H) using the same data in 
(a) and (c), respectively. Here, = 0.4 K and jioHd = 1 T. 



existence of fourfold vertical line nodes. These features can 
be seen more clearly by the contour plot of A^iT, H) in Fig. 
3(b). 

The (f) rotation experiment alone, however, cannot rule 
out the possibility of the horizontal line node as claimed in 
Ref.fTTl In order to solve this issue, we investigated Ce((f>) by 
conic ally rotating H at several fixed 6, where 6 denotes the 
polar angle between H and the c axis. The results obtained 
at 110 mK are shown in Fig. 4(a), where the dashed lines 
are the fit to Eq. ([TJ. In these measurements, the intensity 
of H is adjusted for each 9 to satisfy H/Hc2{6) - 0.2 in or- 
der to avoid unfavorable effects of the tetragonal anisotropy in 
Hc2- The 6 dependence of //c2 was experimentally determined 
from the C^iH) measurements at r=110 mK and is shown 
in Fig. 4(b) by the solid circles. The solid line in the fig- 
ure is a fit to the Ginzburg-Landau anisotropic-effective-mass 
formula for three-dimensional superconductors'^^ Hc2{6) = 
Hc2\\c/{cos^ 6 + {Hc2\\a/Hc2\\c)'^ sin' 5)'^', which reproduces the 
observed 6 variation of //c2 satisfactorily. Figure 4(c) shows 
the dependence of A4 obtained from the results in Fig. 4(a). 
We found that A4 decreases gradually and monotonically with 
decreasing from 90° {H ± c) to 0°. 

Let us compare the experimental results with microscopic 
calculations and discuss the SC gap structure of Celrlns. We 
assume a two-dimensional cylindrical Fermi surface for the 
main Fermi surface of Celrlns with significant /-electron 
contribution.^ The local density of states (LDOS) of the 
quasiparticles (QPs) under the rotating H was calculated by 
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FIG. 4: (Color online) (a) Field-angle (p dependence of C^/T in a 
conically rotating magnetic field with a strength of 0.2H^2(S) at sev- 
eral fixed polar angles 6, measured at 1 10 mK. The dashed lines are 
fits to the data by the expression Co{T) + Ch{T, H)(l -A4 cos 4(p). (b) 
and (c) show the polar angle 9 dependence of //(.2 and A4 at 1 10 mK, 
respectively. The calculated 9 dependence of the zero-energy DOS 
anisotropy [No(H || antinode) - Nq(H || node)] I2Nq(H) is shown in 
(d) and (e) for the order parameters k\ - and k^(kl - kj), respec- 
tively. 



solving the quasi-classical Eilenberger equation self consis- 
tently. Here, we assume the gap function to be A:^. - k^,, and no 
Pauli paramagnetic effect (/i - 0) is considered. The H and T 
variations of the A4 coefficient were then evaluated from the 
LDOS?= and the results are shown in Figs. 3(c) and 3(d) (con- 
tour plot). In the low H (< 0.2//c2) and low T (< 0.12rc) re- 
gion, A4 becomes negative implying that the Cs{(f>) oscillation 
has minima along the nodal directions ((110)). In this regime, 
A4 is approximately proportional to the anisotropy of the zero- 
energy DOS No [IA4I cc 1 -NoiH II node)/No{H \\ antinode)p 
that can be intuitively understood by the Doppler-shift effect 
of the QPs on the circulating supercurrent around the vortices. 
The Doppler shift is given hy 6 = meVp ■ Vs, where nje is the 
electron mass, Vp is the velocity of the QP and Vs is the local 
superfluid velocity that is perpendicular to H. For a supercon- 
ductor with line nodes, A^o is enhanced by the Doppler shift of 
the nodal QPs, giving rise to the //'''' behavior of Ce(//). If H 
is in a nodal direction, then those QPs cannot contribute to A^o 
because 6 vanishes. As a consequence, A^o exhibits an angular 
oscillation with A^o(^ II node) < Nq(H \\ antinode). 
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As T increases, contributions of the finite-energy DOS to 
A4 become relevanti^ia^ The detailed calculations tell us that 
the sign of A4 changes around O-lSTc as shown in Fig. 3(c), 
and Ce(0) takes minima for H \\ antinodal directions in the 
intermediate T range. These theoretical results well repro- 
duce our experimental data as can be seen by comparing the 
contour plots in Figs. 3(b) and 3(d). Especially, the peak po- 
sition and the line of the sign change in A4{T,H) agree re- 
markably with each other The agreement becomes worse if a 
relatively strong Pauli effect is introduced (/U=2 - see Fig. 1 1 
in Ref. l22h . This fact is compatible with the weak Pauli effect 
(yu < 1) expected from Ce{H). Although the low-T Doppler- 
shift-predominant region could not be reached in our experi- 
ment, our data strongly indicate nodes along (110) directions. 

In order to provide further evidence of the d^^i^yi-wawe gap 
in Celrlns, we calculated the 9 dependence of A4. Here we 
compare two types of gap functions k^^ - kf. and (k^^ - k^^k., 
and the results are shown in Figs. 4(d) and 4(e). The lat- 
ter mimics the case of a horizontal line node gap with an az- 
imuthal angular variation of the gap amplitude as claimed in 
Ref. [111. For the c/^.2_,,2-wave gap, a gradual and monotonic de- 
crease of A4 on decreasing 6 from 90° is expected [Fig. 4(d)], 
whereas for the (A:^ - k^?ik^ gap Aa,{0) is predicted to have a 
dip around 9 - 90° [Fig. 4(e)]. The reason why A^iO) shows 
a local minimum at = 90° for the horizontal line node gap 
is because there always exist nodal QPs whose momentum is 
parallel to the field direction, irrespective of (p. In this case, 
the anisotropy in the Doppler shift 5 is strongly reduced, lead- 
ing to the reduction of A4. It is obvious that the k]. - k^ gap 



rather than the horizontal line node gap better reproduces the 
experimental data in Fig. 4 (c). In particular, the absence of 
the dip of A4(0) at = 90° in our results strongly indicates that 
Celrlns does not have a horizontal line node at k~ - 0. These 
observations lead us to conclude that the SC gap in Celrlns is 
of d^2_y2 type. 

In conclusion, we have performed field-angle-resolved 
specific-heat measurements on Celrlns with Tc = 0.4 K by 
conically rotating the magnetic field around the c axis. The 
specific heat exhibits a clear fourfold angular oscillation and 
its dependences on temperature and field is surprisingly well 
reproduced by the microscopic calculations for a c/f2_,,2-wave 
superconductor One of the most important findings in the 
present study is the monotonic variation of the oscillation 
amplitude with the polar angle of a conically rotating field. 
This feature confirms the absence of a horizontal line node on 
the equator and strongly supports the c/i2_,.2-wave gap, as in 
CeCoIns and CeRhlns. The establishment of the identical gap 
symmetry in CeMIns (M - Co, Rh, and Ir) indicates the uni- 
versality of the pairing mechanism in this family and provides 
important hints for resolving the mechanism of the unique su- 
perconductivity. 

This work has been partly supported by Grants-in-Aid for 
Scientific Research on Innovative Areas "Heavy Electrons" 
(20102007, 23102705) from the Ministry of Education, Cul- 
ture, Sports, Science and Technology of Japan. 

Note added in proof. — Recently, Lu et al. also reported 
C(<p) of Celrln5 under pressure down to 0.3 K and the same 
conclusion has been reached.— 
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